Accelerated Cosmological Models in 
First- Order Non-Linear Gravity 

Gianluca ALLEMANDI, 1 '!] Andrzej BOROWIEC, 2 ]] and Mauro FRANC AVIGLIA 1 ' 

^ Dipartimento di Matematica, Universita di Torino 

Via C. Alberto 10, 10123 TORINO (Italy) 
Institute of Theoretical Physics, University of Wroclaw 
PI. Maksa Borna 9, 50-204 WROCLAW (Poland). 
(Dated: February 1, 2008) 

Abstract 

The evidence of the acceleration of universe at present time has lead to investigate modified theories of 
gravity and alternative theories of gravity, which are able to explain acceleration from a theoretical viewpoint 
without the need of introducing dark energy. In this paper we study alternative gravitational theories defined 
by Lagrangians which depend on general functions of the Ricci scalar invariant in minimal interaction with 
matter, in view of their possible cosmological applications. Structural equations for the spacetimes described 
by such theories are solved and the corresponding field equations are investigated in the Palatini formalism, 
which prevents instability problems. Particular examples of these theories are also shown to provide, under 
suitable hypotheses, a coherent theoretical explanation of earlier results concerning the present acceleration 
of the universe and cosmological inflation. We suggest moreover a new possible Lagrangian, depending on 
the inverse of sinh(i?), which gives an explanation to the present acceleration of the universe. 
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I. INTRODUCTION 



Recent astronomical observations have provided strong evidence that we live in an accelerating 
universe. The supernova observation results (see for example 

Q and Q), the observations about 
the anisotropy spectrum of the cosmic microwave background (CMBR) (see for example and 
the results about the power spectrum of large-scale structure [J] have converted cosmologists to 
the acceptance of the evidence of an accelerating universe. 

By itself, acceleration is easy to understand in the context of General Relativity as well as in 
quantum field theory; however problems and doubts on the correct theoretical model to interpret 
observational data arise, owing to the very small but nonzero energy scale which is seemingly 
implied. As a matter of facts, if we believe that we live in a universe which is homogeneous, 
isotropic and accelerating, General Relativity is unambiguous about the need for some sort of 
dark energy source to explain acceleration pj. We are thus faced with some problems concerning: 
(i) the small amount of energy of the vacuum, which is much smaller than we estimate it to be 
(the so-called cosmological constant problem)] (ii) the nature of the dark energy which seems to 
dominate the universe and; (iii) the coincidenceproblem between the actual density of dark energy 



in the universe and the actual matter density 
The real nature of dark energy, which is required by General Relativity in this cosmological 
context, is unknown but it is fairly well accepted that dark energy should behave like a fluid with a 
large negative pressure. This could be explained in vacuum by means of a very small cosmological 
constant, which is in fact related with the energy of the vacuum, or by assuming the presence 
of some matter field, the so-called dark matter or dark energy. The dark energy models with 
effective equation of state w e ff (which determines the relation between pressure p and density of 
matter p) smaller than w e ff < —1 are currently preferrable, owing to the experimental results 
„ t Q. O.fier possibilities include a d_l s CaJar field -*d H , or a phantom 

Q-Q exotic perfect fluids tachvon matter (ic| . 



scalar field 
Born-Infeld Quantum condensate model 



111 ] and a 



a four Fermion model 

One of the first attempts to correctly interpret, from a theoretical point of view, the observed 
data modifying General Relativity (without the introduction of dark energy) was to address the 
cosmological constant problem to somehow allow for the vacuum energy to be large 13]; but this 
has been proven to be not enough to get rid of dark matter and dark energy. 
The other possibility is to assume that we do not yet understand gravity at large scales, which 
means that General Relativity should be modified or replaced by alternative gravitational theories 
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when the curvature of spacetime is small (see for example jl^l). 

Rather than solving the cosmological constant problem or introducing dark matter, we can try to 
explain the current period of acceleration of the universe by a modification of General Relativity 
which leads to modified Friedmann equations (MFR) so that the acceleration kicks in. MFR 
equations of this type arise surely in brane- world models with large spatial extra dimensions [l^ . 
In a completely different framework it has been realized that in the quantization on curved 
spacetimes, when interactions among the quantum fields and the background geometry or the 
self interaction of the gravitational field are considered, the standard Hilbert-Einstein Lagrangian 
has to be suitably modified [153- These corrective terms, which are essential in order to remove 
divergences, are higher-order terms in the curvature invariants such as R, R^R^, R^ ua ^ ] R^ va p, 
RO l R, or non minimally coupled terms between scalar fields and the gravitational field. This 
is a further reason to support the idea that Lagrangians depending on general functions of the 
curvature invariants can provide physically significant models to explain cosmological experimental 
results. It is moreover interesting that such corrective terms to the standard Hilbert-Einstein 
agrangian can be predicted by some time-dependent compactification in string/M-theory (see 



1( 



In particular it has been shown that quantum fluctuations in nearly flat spacetimes 



may induce terms which are proportional to inverse powers of the Ricci scalar invariant for small 
R, while the expansion of the effective action at large curvature predicts terms with a positive 
power of the curvature invariants 3] > [ill • 

As an alternative to extra dimensions, it is also possible to explain the modified Friedmann 



equations by means of a modified theory of four dimensional gravity 



121]. A simple task to 



modify General Relativity, when the curvature is very small, is hence to add to the 



the theory a 



0, 



piece which is proportional to the inverse of the scalar curvature (see 



21] 



n 



agrangian of 



la, ll 



221, 23], 24] and references therein). These theories have been deeply analyzed in 
recent months, in view of their capability of explaining present cosmological acceleration as well 
as early time inflation. As a matter of facts it appears that alternative theories of grav ity provide 



values of w e ft < — 1 and in accordance with the experimental results; see 



and 

We remark moreover that modified theories of gravity, depending on any analytical function of 
the Ricci scalar, have been studied in order to avoid the singularities in cosmological solutions 
j^ . It is worth noticing that black holes in modified gravity are less entropic than in standard 
General Relativity 



IS |. 



It turns out that the metric approach to these Lagrangians leads to modified Friedmann 
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equations which could explain the observed cosmological acceleration of the universe without 
the need of dark energy [lj]. However the metric approach leads to complicated fourth order 



271 ] . It has moreover 



equations that can only be simplified by introducing fictitious scalar fields 
been proven that the aforementioned approach leads to results which are in contrast with the 

n 

solar system experiments |28| and also that the relevant fourth order field equations suffer serious 
instability problems 291. 

It has been shown in [221 ] that, on the contrary, a Palatini variational approach to such modified 
nonlinear gravitational Lagrangians coupled to matter produces second order field equations. 
These equations are equivalent to the standard Einstein field equations only in the vacuum case 
(3^1 or for special (radiating) matter. The modified Friedmann equations obtained in this context 
offer however an alternative explanation for the cosmological acceleration and the approach to the 
de-Sitter spacetime is exponentially fast when the term dominates. These modified Friedmann 
equations are not afflicted by instability problems and they are in acceptable accordance with the 
results of the solar system experiments |2 1| . The Palatini formalism has been further proven to be 
not excluded by electron-electtron scattering experiments provided the physical fields are taken 
into account in a suitable way j^], j^lj . 

On the other hand it is widely believed and accepted that at the very early times the universe 
also underwent an acceleration phase called inflation, the origin of which is still unknown. 
Some authors have proposed that modified gravitational Lagrangians can moreover explain early 
inflation; and |l8|]. The addition of correction terms of the pure-power type R m (m > 0) in 
addition to the j=r term to the standard Hilbert-Einstein Lagrangian, can explain both the early 
time inflation and the current acceleration without instability problems of field equations ^jj] . It 
is there proven that the terms in the Lagrangian with positive powers of the curvature support 
the inflationary period, while terms with negative powers of R are able to explain the current 
cosmic acceleration. 

In this paper we first extend the results previously obtained in (for Lagrangians de- 

pending on general functions of the Ricci scalar invariant in vacuum), to the more general case 
of interaction with matter fields. We consider a first order Palatini formalism for field equations 
of a generalized gravitational Lagrangian in minimal interaction with a matter Lagrangian. The 
gravitational Lagrangian depends on a metric g and a torsionless connection T, assumed to be a 
priori independent. The method generates structural equations for a spacetime with an ensuing 
conformal bi-metric structure. 
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Starting from the request that the connection should be metric with respect to a metric h (suitably 
defined by dynamics) we are able to obtain generalized Einstein field equations. The metric h 
turns out to be related to the metric g by a conformal transformation. These general results 



extend those of the vacuum case 



301 ] : the universality of the Einstein field equations does not 



hold anylonger, apart from the case of a purely radiation universe (in this case the stress-energy 
tensor is traceless and the structural equations for spacetime are the same as in the vacuum 
case). The universality property was investigated not only at the level of equations of motion, 
but also at the level of the definition of energy and conservation laws for the gravitational field. 
This problem was in fact analyzed in the paper 32], where it was shown that universality holds 
also for the gravitational energy-momentum complex, which turns out to be essentially the same 
as in the linear case. The formalism, developed in this paper to obtain the generalized Einstein 
field equations, is very general and in forthcoming papers (under completion) it will be moreover 
generalized to theories depending on higher-order terms in the Ricci-squared curvature invariant 
[33J, which have been already analyzed in the vacuum case in [3J|. 

For cosmological applications, we shall substitute the general Robertson- Walker metric g in the 
generalized Einstein field equations, assuming the stress-energy tensor to be a perfect fluid tensor. 
Taking into account the conformal relation between the (physical) metric g and the conformal 
metric h, we obtain generalized Einstein equations which reproduce a modified Hubble constant, 
depending on the conformal factor. We stress that, owing to the conformal transformation, also 
the metric h results to be Robertson- Walker. The formalism is fully and analytically applicable 
under the request that the structural equation can be solved so that both f(R) and f'(R) can be 
expressed in terms of the trace r = Tn V g^ v (see also 3]). Otherwise, numerical and approximation 
methods are always applicable. These alternative theories are very rich in their structure and they 
can provide models of multi-universe cosmologies and can be suitably modified to fit experimental 
data. 

Examples can be provided by choosing particular expressions of f{R). We treat here the simplest 
cases of f(R) = (3R n (pure power) and f(R) = a ln((3R) (logarithmic Lagrangian), which are 
very significant in view of their cosmological applications. The modified Hubble constant and the 
first order field equations are obtained, both for h and g. We propose moreover a new possible 
Lagrangian, which deserves interesting cosmological applications, containing terms in the inverse 
of the hyperbolic sinus of the Ricci scalar. 

The example f(R) = j3R n is remarkably important. As a matter of facts each analytical La- 
grangian in R can be locally approximated by means of a Taylor polynomial expansion and further 
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on polynomial Lagrangians in R can be (under suitably hypotheses) locally approximated by 
f(R) — f3R n , where n will be negative for small R and positive for large R. In the aforementioned 
case of f(R) = (3R n we specify the stress-energy tensor for the cases of dust, radiation and 
vacuum universes. In the dust universe and for the particular choice of zero space curvature in the 
Robertson- Walker metric we can solve exactly the differential equations and we obtain the value 
of the acceleration parameter depending on the exponent n (for particular values of n). It results 
that acceleration is supported in the cases n < and n > |. These models seem to be able to 
explain both the current universe acceleration and early time inflation, taking into account leading 
terms of polynomial-like Lagrangians in each epoch of the universe. Polynomial Lagrangians are 
also discussed in order to treat radiation universes. It results that present time acceleration can 
be predicted and a critical value of the scale factor appears, which rules the transition from an 
accelerating to a decelerating epoch of the universe. 

Logarithmic Lagrangians in the Ricci scalar give an interesting, exactly solvable, model also for 
the case of pure radiation matter, supporting the present time acceleration of universe (see also 
Q and Q). 

Lagrangians containing terms proportional to the inverse of the hyperbolic sinus of the Ricci 
scalar provide very complicated modified Friedmann equation which cannot be solved analytically. 
However these equations could give an explanation for the current acceleration of the universe in 
the appropriate limit and they can be useful for further numerical analysis or finer approximation 
techniques. 



II. f(R) GRAVITY 

We begin with considering on a 4-dimensional Lorentzian manifold (M, g) the action 

A = 4 rav + A mat = J (^gf(R) + 2KL mat )d 4 x (1) 

where R = R(g,T) = g a @ R a p(T) and R^iT) is the Ricci tensor of any independent torsionless 
connection T. The gravitational part of the Lagrangian is controlled by a given real analytic 
function of one real variable f(R), while ^fg denotes the scalar density | det || g^ u |||2 of weight 
1. The total Lagrangian contains also a first order matter part L mat functionally depending on 
matter fields together with their first derivatives, equipped with a gravitational coupling constant 
k = 8ttG. 
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Equations of motion ensuing from the first order a la Palatini formalism are (see 



22, a3 



f(R)R iflv) (T)-±f(R) gixu = kT„ u (2) 
VK^gf'W) = (3) 

where T^ v = — %= ^ mat denotes the matter source stress-energy tensor and V r means covariant 
derivative with respect to T. In this paper the metric g and its inverse are used for lowering and 
raising indices. 

We shall use the standard notation denoting by Ruv) the symmetric part of R^ u , i.e. R(^u) = 
\{R^u + Rv/j,)- In order to get © one has to additionally assume that L mat is functionally inde- 
pendent of T; however it may contain metric covariant derivatives V 9 of fields. This entails that 
the matter stress-energy tensor T^ u = T fiu {g^ ,V 9 ^) depends on the metric g, the matter fields 
denoted here by together with their metric covariant derivatives. Physically speaking we are 
assuming that matter fields are minimally coupled to the gravitational field g. 
From Q one sees that y/gf'{R)g^ u is a symmetric twice contravariant tensor density of weight 1, 
so that if it is not degenerate one can use it to define a metric hav such that the following holds 
true 

y/gf{R)gr = VhhT (4) 

This means that the two metrics h and g are conformally equivalent. The corresponding conformal 
factor can be easily found (up to a possible sign) to be f'(R) (in dimM = 4) and the conformal 
transformation results to be: 

h/ju, = f{R)g l xv (5) 

Therefore, as it is well known, equation (j3J) implies that V = Ti,c(h), i.e. the Levi-Civita connection 
of h, and ^(^(r) = R^ u (h) = R^. We should emphasize, however, that scalar R = g^R^ih) 
is not a Ricci scalar of the metric h. In what follows we shall call it a generalized Ricci scalar and 
sometime for shortcut simply a Ricci scalar. 

Equation can be supplemented by the scalar-valued equation obtained by taking the g-trace of 
(J2J), where we set r = trT = g^T^: 

f'(R)R-2f(R) = KT (6) 

We obtain that equation © controls solutions of equation (|2~]l. We shall refer to this scalar- 
valued equation as the structural equation of the spacetime. More precisely, for any real solution 



R = F(t) of © we have that f{R) = f(F(r)) and f'(R) = f (F(r)) can be seen as functions of r. 
For notational convenience we shall use the abuse of notation /(r) = f(F(r)) and /'(r) = f'(F(r)). 
Now we are in position to introduce the generalized Einstein equation under the form 



For the matter-free case 



with h^ u defined by © for a given g^ u and (see also 
we find that R = F(0) becomes a constant, which implies that the two metrics are homotetic; this 
property further implies that equation is nothing but Einstein equation for the metric g and 
it is almost independent on the choice of the function f(R). This is nothing but the universality 



property observed in 30] . Also Einstein equation with cosmological constant A can be recasted into 
the form © if we chose f(R) = R — A. These properties justify the name of generalized Einstein 
equation for I©. In the presence of matter ((TJ) expresses instead a deviation for the metric g to be 
Einstein, apart from the case when r = where universality still holds true 3oT ]. 
It should be also noticed that equation @ has, in general, many real solutions, especially when 
f(R) is a polynomial function of higher degree. Therefore, /(ii)-gravity supports the idea of multi- 
universe interpretation. Moreover, the number of possible universes is dynamical, since r turns out 
to be a function of the FRW scale factor a(t) (see below) except for a radiation dominated period. 
Another special property one wants to emphasize is that we can make f'(R)R — 2f(R) to be any 
function we wish. This is, of course, due to the fact that the non-homogenous linear ODE 

f{R)R-2f{R) = 4>{R) (8) 

has a solution in the form f(R) = R 2 J ds^^-. The choice of f(R) allows to design a specific 
cosmological scenario and to adjust the model to fit concrete experimental data. 

III. FRW COSMOLOGY IN f(R) GRAVITY 

For the cosmological applications one has to choose the metric g to be the Friedmann-Robertson- 
Walker metric, which (in spherical coordinates) takes the standard form: 

g = -dt 2 + a 2 {t) - _^ Kr2 dr 2 + r 2 (d9 2 + sin 2 (0) dy> 2 ) (9) 

where a(t) is the so-called scale factor and K is the space curvature {K = 0, 1, —1). Another main 
ingredient of the cosmological model is the perfect fluid stress-energy tensor 

T^ u = (p + P)u^u v + pg^ v 



1 /IV 



p o 

u T=Kr* 



where p is the pressure, p is the density of matter and is a co-moving fluid vector, which in a 
co-moving frame (u M = (1,0,0,0)) becomes simply: 

In n o o\ 


pa 2 (t)r 2 
\0 pa 2 (t)r 2 sin 2 (<9) J 

Later on we shall assume the standard relations between the pressure p, the matter density p and 
the expansion factor a(t), namely: 



(10) 



p = wp , p = rja 



-3(1+™) 



(11) 



where particular values of the parameter w G {—1,0, ^} will correspond to the vacuum, dust or 
radiation dominated universe. Exotic matters, which are up to now under investigation as possible 
models for dark energy, admit instead values of w < — 1. These expressions (|11|) follow from the 
conservation law of the energy-momentum V^T^ = and consequently the continuity equation 
should hold: 



p + 3H(p + p) = 



(12) 



where H = - is the Hubble constant. The requirements © and @ fix h to be conformal to g and 
in fact equal to: 

1 



h 



f'(r){-dt 2 + a 2 (t) 



; dr 2 + r 2 (d9 2 + sin 2 {9)d<p 2 ^J } 



(13) 



.1 -Kr 2 

where r = T^ u g^ = 3p — p is a function of time through its dependence on the scale factor a(t): 

t=(3w- l)r][a(t)]-^ 1+w ^ . 

Substituting all necessary ingredients into the generalized Einstein equation we obtain the 
following 



^ (fc) _=/^__» [a § + >^_(»)* 



(14) 



for the 00 component, while for the 11 component we have: 
T-iTr 2 l f(r) + 2K P 1 



Ru(h) 



2b 



S b a b A /d\2 K 
a b a b \aJ a z 



(15) 
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Here for simplicity we have introduced the shortcut b(t) = /'( T )- Of course for the usual Einstenian 
case one has b(t) = 1. Combining the last two equations we can obtain an analogue of the Friedmann 
equation under the form 



a b \ 2 K k r f(r) + kt 



a + 2b) + a 2 36 



P + 



2k 

b 



(16) 



which can be seen as a generalized definition of a modified Hubble constant H = ( ^ + 55 J , which 
takes into account the presence of the conformal factor 6, entering the definition of the conformal 
metric h, which in turn generalizes the results of [21]]. 

We stress that, owing to the conformal transformation between h and g, the generalized Ricci 
scalar can be generally expressed from (|14|) and (|15|) as follows: 

R = = 2/(r) fe + KT = R(h)b (17) 

which reproduces General Relativity in the particular case f(R) = R, as we shall see later on, but 
provides relevant deviation from General Relativity in the other cases. We remark that R(h) is the 
true Ricci scalar of the metric h. As it has been already remarked in the metric h itself is a 
Robertson- Walker metric (see p3[l. 

h = e(-dP + A 2 {t) - —^-3 dr 2 + r 2 [dd 2 + sin 2 (g)dy 2 ) ^ (18) 

with a new cosmic time dt = y^\b\dt and a new scale factor A = ^/\b\a, where e = ±1 corresponds 

to positive or negative values of b. The generalized Einstein equation (J7J) can be also calculated in 

(t,x t )- coordinates. This is equivalent to the assumption that the metric h is a physical one (i.e., 

that we can use conformal frame instead of the original Einstein frame). We obtain in this case: 

-m + 2n P _ A 
e 26 2 " A [ 9J 

for the 00 component while for the 11 component we find 

where A denotes now the differentiation with respect to the new cosmic time t. An additional 
factor I on the l.h.s. appears due to the fact that dt 2 = ^dt 2 and a 2 = j^^ 2 - Now the analogue 
of the Friedmann equation takes the form 

~ 2 K /(r)+«(p + 3p) 

H = ~A^ + e 6P (21) 

with H = 4 being the Hubble constant of the conformal metric h. This expression is much simpler 
and recalls the standard Friedmann equations for the Robertson- Walker metric. 

10 



A. Examples: f(R) — (3R n and polynomial Lagrangians 

We consider, as a class of particular examples, the class of linear Lagrangians in an arbitrary power 
of the curvature invariant R. The importance of such models can be considered and understood 
in connection with the main effort of modifying General Relativity by means of alternative La- 
grangians, which are able to explain the experimental data in some limit. Up to now (and up to 



our knowledge) polynomial-like and logarithmic- 
invariant have been considered (see for example 



ike Lagr ang ians in the generalized Ricci scalar 

3, Q, 0, Q, Q, Q, Q 



and refer- 



ences therein). These particular Lagrangians recover General Relativity in some approximation, 
which physically speaking means at some age of the universe, and reproduce modified Friedmann 
equations in some other limit. These modified Friedmann equations are able to give a possible 
theoretical explaination to the experimental results (such as inflation and present acceleration of 
the universe). 

If we narrow down our researches to the case of polynomial Lagrangians, each term of the La- 
grangian behaves as a leading term at some particular age of the universe and it is able to reproduce 
at a convenient order of approximation the experimental data. Moreover any analytical function 
can be approximated by means of its Taylor polynomial expansions in the limits of physical rele- 
vance and it consequently behaves like a polynomial-like Lagrangian. 

In this framework it is hence worth analyzing exactly and analytically Lagrangians of the type 
L(R) = f3R n ^/g (with arbitrary possibly non-integer n) representing single terms of a more phys- 
ical and general polynomial or polynomial-like expression for the Lagrangian. Integer values of 
n assume a fundamental role in this context representing, for any analytical function, the terms 
deriving from a Taylor expansion. The advantage which derives from this particular class of La- 
grangians is that they are easily and exactly solvable and they provide coherent models for the 
universe acceleration and for the early time inflation. 

We stress once more that we are specializing to the case of four dimensional spacetimes so that 
9^ v 9\w — 4- We assume then a Lagrangian of the form: 

f{R)y/g = (iR n V9 W + 0; n e R; n ± 0, 2) 
(see also Q] and 0i). The algebraic field equations (J2J) are consequently the following: 



(3R n 1 nR^y — -Rg^ 
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while the structural equation © becomes: 

R" 



KT 



P(n 



(22) 



This expression implies that i lie case n 1 is singular (see e.g. ^|). Taking into account the 
expression (jlOj) for the stress energy tensor and the dominant energy condition we should impose 



r = T^ v < 



(23) 



Equation (|22|) forces R n to be positive definite for even integer values of n and we should fix the 
dimensional constant (3 (the dimension of is the same as the dimension of R l ~ n ) in front of the 
Lagrangian to be : 

j > for n < 2 
\ < for n > 2 

Let us introduce, for convenience, a new re-scaled dimensional coefficient (3 = 0(2 — n) which is 
positive for even n and arbitrary otherwise. In the last case, i.e. for n odd, one has signi? = sign/3, 
where R = (— Kr/0)n. Instead, for n even one has two different real solutions of ()22j) . namely: 



R± 



- KT \ n 



ft 



(24) 



This implies that in any case R will be proportional to a well-defined power of a, which is exactly: 

3Q+1) 

R ~ a n This expression shows that only in the vacuum case the model approaches a de- 
Sitter (anti de-Sitter) universe. 

The requirements (JHJ) and (@J fix h to be conformal to g and in fact equal to: 

ne 



h 



n 



-f3n (—kt)~ 



(25) 



where we remark again that e = signi? = 1 for odd values of n and, on the contrary, e = ±1 in 
accordance with the choice of the solution in (|24|) for even values of n. We are now able to calculate 
the modified Friedmann equations and the modified Hubble constant from H16|) : or, in a different 
but completely equivalent way, by inserting (|25ft into 0. From the structural equation Q we can 
calculate: 

r m = 

and we obtain that the Hubble constant for the metric g can be locally calculated to be: 



H 



2ne 



3(3u; - l)[3w(n - 1) + (n - 3)] 



-KT 







K 



2n 



3w(n — 1) + (n — 3) 



(26) 
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The deceleration parameter can be obtained from the Hubble constant by means of the following 
relation: 

*>--(-j§)--UiU 

and from (|26jl it results, in the case K = 0, to be formally equal to: 

s 3(l + w)-2n , , 

q(n,w) = -± ^ (28) 

We say formula (|28|) to be formal since, at the moment, we do not know about its effective solution 
realizations. This will be the subject of our investigations below. 

The effective w e tf can be obtained (as in by means of simple calculations from (|28|l and it 
results to be, for this theory: 

2 , . 1 1 w 

Weff = - q (n, W )- 3 =-! + - + - (29) 



We remark that the range of —1.45 < w e ff < —0.74 for dark matter, stated in [3j, can be easily 
recovered in our theory by choosing suitable and admissible values of n. 

If we consider expression ()28|) for the deceleration parameter q(n,w) we see that, from the 
definition of w e ff in (|29jl. an accelerating behaviour of the cosmological model requires that the 

effective w e jf < — |. This bound to the value of w e ff can be alternatively seen in terms of w. 

2n 
3 



This states in particular that to obtain accelerated universes we should else impose w < Tp — 1 for 



n > or w > y - 1 for n < 0. 

For the case n = 1 this imposes an upper bound for the acceleration to the value of w < w cr a = — 1, 
reproducing the well-known results of General Relativity; choosing alternative Lagrangians with 
n > implies that this limit is shifted to the value w cr u = — 1. In the case n < we have 
no longer an upper bound, but a lower bound for w > w cr u = ^p — 1; dust and radiation matter 
satisfy that condition by definition. 



Consider now, in more detail, the case of pressure-free (dust) universe, i.e. p = 0, so that r = 
—pit) = ~aT(jj ( we remark once more that 77 > 0). The resulting generalized Einstein equations 
(|26jl which derive from (j2J produce modified Friedmann equations in the matter universe case 
considered here. The resulting expression for the Hubble constant of g can be calculated from (|26|) 
to be (for 



a 2 (t) 
a 2 it) 



eniijn) 



6Kn 2 



(3-n)[/3a 3 ]^ a 2 (t)(n - 3) 2 



(30) 
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which reproduces, as it should be expected, the standard Friedmann equation in the very particular 
case n = 1, and j3 = (3 = 1. We remark that the above expression (j3()j) represents the square value 
of the Hubble constant and it becomes singular at n = 3 in the case of dust universe, under analysis. 
We restrict now ourselves to the case of K=0. We have to require that the expression on the right 
hand side of (|3*U|) should be positive. 

We start discussing integer values of n, which are very relevant for our analysis, as we stated before. 
In the case of odd integer values of n (e = 1) we should require (3 > for 0<n<2or2<n<3, 
i.e. essentially in the case n = 1 of standard General Relativity. When n is still odd but (3 < one 
gets re < or n > 3 instead. 

In the case of even values of re, we see that (3 is always positive. Therefore, as we required before, 
we should fix e = — 1 for the following values of n: re < or n > 3. Thus the solution R + of ()24|) 
gives no relevant contribution to the solutions of (|30|) . owing the singularity in the case n = 2. 
Finally, one should notice that (3 > for n G (0, 2) U (3, +oo) and negative otherwise, i.e. for 
n E (— oo, 0) U (2, 3) independently of the parity of n. 

We remark however that we do not need to assume that n is an integer and, in fact, it can be 
a priori any real (or rational) number. In any case, however, we should require equation ()22j) to 
be definite and morever we should impose the positivity of (|26|) , Suitable values of (3 have to be 
chosen and the analysis can be carryed over, following the headlines of the analysis previously done 
in the case of integer values of n. Specializing again to the case of dust universe with K = 0, in 
the case of real, non-rational n one has to assume (3 > and this imposes n G (0, 3). 
The modified Friedmann equations, for dust universe in the case K = 0, can be integrated and we 
easily obtain (apart from integration constants, which can be forgotten in our analysis): 



a(t) 



3e 


n 
3 


KT] 


2re(3 - n)_ 




IP] 



tT (31) 



Now the deceleration parameter can be obtained again from (|28|) : we find: 

. , 3 — 2n . , 

(Kt) = — (32) 

for the particular solutions (|3*T|) corresponding to the case K = 0, selected by (|3*U|) . This implies 
that we obtain accelerated solutions in the case: 

3 

q(t) <0^n<0 orn>- (33) 

corresponding to alternative theories of gravity with Lagrangians depending on inverse powers of 
the generalized Ricci scalar or on terms with powers higher than |, owing to the restriction imposed 
for the positivity of ()3U|) . 
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If we consider the case of pure-radiation universe, corresponding to r = 3p — p = with 
p = p(t) = -Jrijji w e have that the equation (JJJ) is undefined, since (j2*2*|) implies R n = and 
/'(t) = 0. This entails that our formalism fails to cover this case which has to be treated 
differently (see below). Therefore, formula (|28|) for w = | is purely formal. 



In vacuum dominated universe, corresponding to the case w = —1, we obtain from l)26|) that the 
Hubble constant is: 



H 2 _ a 2 (t) 



2,n - r 



P 



(34) 



o 2 (t) 



a 2 (t) 12 

which results to be positive and well-defined for any n (n / 2) with a suitable choice for (5 and e. 
It follows that in the case of a vacuum universe, and for K = we have that: 



a(t) = eV «L /a J (35) 

corresponding to a de Sitter universe (with a power-law expansion); the deceleration parameter is 
in this case q(t) = — 1 supporting cosmological acceleration. 

A similar analysis can be performed under the hypothesis that the metric h is the physical metric 
and it is spatially flat {K = 0). The deceleration parameter for the metric h can be formally 
calculated to be: 

-/(t) + 2k P 2n-3(l + w) 2n-3(l + w) 

f(T) + K(p + 3p) n-3 + 3(n-l)w n(l + 3w) - 3(1 + w) 1 j 



which gives acceleration q(—l,n) = —1 for any value of n in vacuum dominated universes and 

2ra-3 
n— 3 



c/(0, rt) = 2 "_o 3 < in dust dominated universes provided that | < n < 3 (n 7^ 2). Finally, for pure 



radiation we have q(^,n) = 1, i.e. one gets deceleration for any value of n. 
The effective value of w can be also calculated in this case from (|29[) : we get: 

4n-9(l + w)+ra(l + 3tt;) 
Weff ~ 3[n(l + 3w)-3(l+w)} [ ' 

which allows us to make a comparison with the experimental data. 

These results are alternative to the case described by formula (|28|) when g is the physical metric. 
As usual in alternative non-linear theories, we do not know a priori which is the physical metric: 
discussions about the physical interpretation of q and h both from the mathematical and physical 

i. g . 0, 0h 



viewpoint are up to now open; see, e.g. [27| . |3J| and [37 1. 
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1. Polynomial Lagrangians in the generalized Ricci scalar 



As we told before theories with power Lagrangians in the generalized Ricci scalar can be considered 



as approximations of more physical polynomial-like Lagrangians 



3 



of the type: 



f(R) = R+ - ^— - + — R m (38) 

(here both n > and m > 0, with m/2 and n ^ —2). 

We consider the case of dust universe w = 0. In the limit of small or large curvatures, correspond- 
ing to the cases of present time universe and early time universe, we obtain from the structural 
equations that the leading terms are respectively: 

r r^q^^ = kt 

\ R — > co =^ -/3i? m = KT 
From (|32|) we deduce that polynomial Lagrangians provide an explanation for early time inflation 
assuming that m > | and they can provide an explanation to present time cosmic acceleration 
assuming that some inverse power of the generalized Ricci scalar is present in the Lagrangian (i.e. 
a / 0). 



This result reproduces previous resu 



in the metric formalism (in 



13] and 



ts which have been obtained in a different framework both 



19]] 



and in the Palatini formalism ( [211 ] and |22|). We stress 



3- 



moreover that terms like (3R m (with m > 0) are related to the so-called Starobinsky inflation 



2. Polynomial Lagrangians for radiation universes 

We can otherwise consider, for the case of radiation universes w = 3, a Lagrangian of the form: 1 

f(R) = aR+ -^—R m (where w = -; r = 0) (39) 
m — 2 3 

which in 

Q and has been already examined in the case m = 2,t / 0. It has been there 
proven that, in the metric formalism such Lagrangians support inflation, while in the Palatini 
formalism they provide explanation for the present time acceleration [2l| . 

The structural equation Q admits besides the trivial solution R = also a non-trivial solution 
1 

R = !- m , where e = 1 for m odd and e = ±1 for even values of m. This correspond to de-Sitter 
and anti-de Sitter universe for positive and negative values of R respectively. 



m # 1,2 
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In the case of the obvious solution R = we have that f(R) = 0, f'(R) = a, and consequently: 

3a 3a 

The deceleration parameter can be formally calculated to be: 

Q(t) = m Y\ a -2 ( 41 ) 

In the particular case K = we have a = ^^\/2t and the deceleration parameter can be easly 

calculated to be q(t) = 1, describing a decelerating universe for this solution (R = 0). 

If we shift to the non trivial solution R = e (^J m 1 we now see that the structural equation gives: 

f{R) = e^-^(ir X (42) 



m-2 \p 

so that the Hubble constant for the metric g can be calculated to be: 

i 

2 Krjim-2) _ 4 m-2 + a fa\ m - 1 _ 2 „ T . _ 2 a -4 /aa\ 

R = — -a + e ; r - - Ka A = S - ^ + Aa 4 (44) 

6a(m-l) 12a(m-l)V/3/ 

where for notational convenience we have introduced the parameters A = g^^Zi) an d ^ = 

m-2+q ( a \ m ~ 1 



e 12 ^(m"-i) (^f J ' wn i cn we require to be positive. 
The deceleration parameter can be obtained by means of formula (|27|): 

Aa- 4 -S A-Sa 4 
m) = + = A - Ka 2 + Sa 4 (45) 

In the limit of large density matter p, corresponding to early time universes, we obtain q = 1; 
on the contrary at late times we have q = — 1 when p results to be very small and forgettable. 
This corresponds to a presently accelerating universe. We remark that the critical value q = 
(corresponding to a change from a deceleration epoch to an accelerated epoch) is met when the 
radius attains the critical value: 



i 



A i 2enr] (m - 2) ( (3\ ™-! 



y m — 2 + a \a / 

We remark that both asymptotic values of q and the critical value a c do not depend on the value 
of the spatial curvature K of the spacetime under consideration. 

This simple example illustrates two important properties of polynomial Lagrangians, which are full 
in physical significance: 
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- there may exist two (or more) parallel universes, corresponding to different solutions of the 
structural equations of the same Lagrangian and matter source. These different solutions 
provide models for different cosmologies; 

- there exist solutions which provide a smooth transition from deceleration epochs to acceler- 
ated universes. This happens in correpondence with some critical value of the cosmic radius 
a c {t) and it is in relation with the so-called cosmic speed-up [~Lt1 |. 



B. Another example: f(R) — a\n(R) Lagrangians 



We chose as a further particular example the case f(R) = aln((3R), which is relevant since log 
arithmic terms in the Ricci scalar are induced by quantum effects in curved spacetimes; see |l£ 
and Notice that the corresponding dimensions should be [a] = t4t = [R]. For notational 

convenience we will fix (3 = 1 in units such that dimensions remain correct. 
We obtain from the structural equation © that: 

f Sir) = 9=*z 



{ b(t) = f{r) = ae— 

Substituting this expressions into Q16|) and performing straightforward calculations the explicit 
expression for the Hubble constant can be easly obtained: 

1 



H 



1 



3k(1+w) , 



la 



3(w + 1)kt + a(3w - 1) K 

~~1 



12a{3w - l)e—° 



(46) 



where we recall that r can be expressed in terms of a(t) as: 

T = 3p - p = r](3w - l)[a(i)]- 3 ( 1+ ^ 



(47) 



We remark that the Hubble constant can be specialized to the case of dust, radiation or vacuum 
universes by a suitable choice of w (w = — 1, 0, 3). It follows from formula (|4*6]> that the pure 
radiation case is singular and cannot be treated in our formalism. 
In the case of dust universe (w = 0) we have that: 



1 



a + 3ko <x+n P 



K 
72 



(48) 



12a 

We consider the limits of small and large energy density in the universe, which are respectively 
a » up and a << up. In the case of small energy density, corresponding to late times of the 
universe, we have that: 

r2 _Ve 



12 



(49) 



18 



which reproduces a power-law expansion with a universe approaching a de-Sitter universe. On the 
other hand, in the case of large density matter we have that, apart from constant and positive 2 
factors, the following holds H 2 ~ e p and this surely supports decelerating cosmological models. 
The deceleration parameter can be moreover calculated formally. It turns out to be, from l|27|): 



q(t,w) 



9k(1 + w) 2 t 

■1 H h 

3k(1 + w)t — 4a 



(50) 



3k(w + 1)t[3(w + 1)kt + a(-3w - 7)] 2K 



24a 2 (3w - l)e^ 



+ 



3(w + 1)kt + a(3w - 1) K 
I2a{3w - l)e*%r a 2 



which can be suitably approximated in the cases of small and large energy densities. 
We obtain that in the limit of large energy density, corresponding to early time universes, we will 
have q(t) ~ 3(u> + l) — 2, which provides decelerating universes, apart for the case of vacuum space. 
The limit of late time universe, namely for small energy density r, provides as expected from (|48fl 
a de-Sitter like universe with q(t) = —1, which always provides an accelerating cosmological model. 
The analysis performed above implies that, in the Palatini formalism, Lagrangians proportional to 
the logarithm of the Ricci scalar provide cosmological models without an inflationary epoch, while 
these models are able to explain the current acceleration of the universe; see also |3J3] and 3j|. 



C. A possible new cosmological model: f(R) = R — sin ^" fl ) 

We consider moreover an alternative Lagrangian containing a term proportional to the inverse of 
the hyperbolic sinus of the generalized Ricci scalar: 

f(R) = R- . *** (51) 
v ; smh(i?) v ; 

where coefficients have been chosen to simplify future equations and reproduce General Relativity 
in the appropriate limit. Structural equations © hold in the form: 



R cos h(.R) + 2 si nh(.R) 
sinh 2 (i?) 



R-Qa "^ 11 ^,2 „ =~kt (52) 



We consider the limit of the above structural equation for the case of small values of R (see also 



3 



2J and 



221]), Solutions can be found under the form: 



-kt ± v^kt) 2 + 72a(l + a) 

R = WTa) (53) 



the positivity is necessary and it follows from the positivity of H 2 



19 



To reproduce, in the limit | kt |>> a (see 



the result for General Relativity R 



l+a 



KT 



we have to choose the plus sign in 1)53)) . We remark that in the case a = General Relativity is 
exactly recovered, as it should be expected. This choice implies that in the limit of late universe 
a >>| kt I we will have: 



R 



18a 

[l + a) 



so that deviations from Einstein gravitational theory are large and the universe approaches in that 
limit a de-Sitter universe, dominated by a power-law expansion. 

We obtain, in this case, that equation 1)16)1 for the Hubble constant cannot be solved analytically; 
neither the structural equation 1)53)1 provides an explicit analytical expression of R in terms of r. 
We can however write the system in an implicit form: 



-9(w + l)ansmh 2 (R)(l + cosh 2 (fl)) 



(smh^-R) + 6acosh(i?))[sinh 3 (i?) + 6a(R + Rcosh 2 (R) + sinh(i?) cosh(i?))] 



Ksinh 2 (R) 



R 



6a 



3[sinh 2 (ii) + 6acosh(i?)] 



P + 



sinh(R) 



+ KT 



2k 



„ Rcosh{R) +2sinh(i?) 

R — 6a ^ = —kt 

sinh 2 (ii) 



(54) 



(55) 



where the solution R = F(t) of 1)55(1 . substituted into 1)54)1 . gives a formal solution to the system. 
This expression could be possibly useful for some further numerical analysis. 

Considering the leading terms in the limit of small R, corresponding to the present universe we 
are interested in, the above system can be solved in a rough approximation. We obtain from the 
structural equation 1)51-1)1 that R = and substituting in the approximated expression for H 2 we 
obtain: 



H 2 ~ [3w + 4] 5 



3a \ 

KT J 



a > 



(56) 



so that the current acceleration of the universe is explained by means of this model in corre- 
spondence with small values of R. As a matter of facts, apart from constant factors, we have 
that H 2 ~ t -1 and consequently for the most interesting case of dust universe we have that 
H 2 ~ a 3 (t), describing an accelerating universe. In this cosmological model we will have that 
the effective w e ff & t present time, which means in the limit 1)56)) . can be simply obtained to be 
w e ff = — 2 — w, which goes in the direction of the experimental results of [3j. 
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IV. CONCLUSIONS AND PERSPECTIVES 



In this paper we have analyzed alternative theories of gravity, the Lagrangian of which is a 
general function of the generalized Ricci scalar R, constructed out of a dynamical metric g and 
a dynamical connection V. The Palatini formalism provides first order field equations for the 
metric and the connection T. A structural metric h is introduced, such that the connection 
results to be the Levi-Civita connection of h and h is consequently conformal to g. Spacetime 
geometry is thus defined by means of generalized Einstein equations and it is in fact mutuated 
from structural equations which reproduce the standard Einstein spacetime with suitable choices 
of the parameters. 

To treat explicitly cosmological models we choose g to be a Robertson- Walker metric. This allows 
to obtain modified Friedmann field equations and a modified Hubble constant related to the 
conformal transformation between g and h. The metric h results to be FRW, too, so that it can 
be conveniently considered as a physical metric in place of the original g. 

If we moreover specialize to the pure-power case f(R) = f3R n (with n an arbitrary real exponent) 
we have seen that, with suitable choices of the parameters involved, these models are able to 
explain the current acceleration of the universe. We obtain that polynomial Lagrangians in the 
generalized Ricci scalar provide an explanation for both present acceleration and inflation of the 
universe in suitable limits 

Considering instead the case of logarithmic Lagrangians in R we have that in the limit of small 
matter density, corresponding to late time universe, the solution approaches a de-Sitter universe 
with exponential power law expansion. This results into a possible model for accelerating universe. 
We also propose a new possible cosmological model inlolving the inverse power of the hyperbolic 
sinus. Even if it is not analytically solvable it however provides, in a suitable limit, an explanation 
for the current acceleration of the universe as it should be expected. 

The formalism developed here is very general and mathematically well defined. It provides useful 
physical results in the specific cases we have considered. Of course, a fully covariant satisfactory 
Lagrangian should be a more complicated polynomial-like expression involving more than one 
power R m and more than one inverse power R~ n , besides a (possible) logarithmic factor ln(/3R) 
and terms proportional to the inverse power of sinh(i?). In this case, however, the non-linearity 
of the structure equation does not allow in general a simple analytical resolution as we did in the 
specific examples, so that numerical or approximation techniques have to be invoked for generic 
Lagrangians. 
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Moreover, we stress that further applications of the formalism developed here can be ex- 
tended to cover also generalized Ricci-squared theories, which will be analyzed in a forthcoming 
paper |33|. 



V. ACKNOWLEDGEMENTS 

We are very grateful to Prof. S.D. Odintsov for useful discussions and very important suggestions, 
concerning the physical properties of the cases considered, and for suggesting us some of the ex- 
amples. 

This work is partially supported by GNFM-INdAM research project "Metodi geometrici in mecca- 
nica classica, teoria dei campi e termodinamica" and by MIUR: PRIN 2003 on "Conservation laws 
and thermodynamics in continuum mechanics and field theories" . Gianluca Allemandi is supported 
by the I.N.d.A.M. grant: "Assegno di collaborazione ad attivita di ricerca a. a. 2002-2003". 



[1] Perlmutter, S. et al. 1999, ApJ, 517, 565 (astro-ph/9812133 ); Perlmutter el al. Nature 404 (2000) 955 

Astroph. J. 517 (1999) 565. 
[2] Riess, A. G. et al. 1998, AJ, 116, 1009 (astro-ph/9805201). 

[3] C.L. Bennett et al., Astrophys.J.Suppl. 148 (2003) 1 (astro-ph/0302207); C.B. Netterfield et al., As- 
trophys.J. 571 (2002) 604-614 (astro-ph/0104460); N.W. Halverson et al., Ap.J. 568, 38 (2002), (astro- 
ph/0104489); A.H. Jaffe et al., Phys. Rev. Lett. 86, 3475 (2001), (astro-ph/0007333); A.E. Lange et 
al., Phys. Rev. D 63, 042001 (2001); A. Melchiorri et al., Ap.J. 563, L63 (2000), (astro-ph/9911445); 
Spergel, D. N. ct al. 2003, ApJS, 148, 175. 

[4] Verde, L. et al. 2002, MNRAS, 335, 432. 

[5] S. M. Carroll, Measuring and Modeling the Universe, Carnegie Observatories Astrophysics Series Vol. 
2, ed. W. L. Freedman (astro-ph/0310342). 

[6] R. R. Caldwell, R. Dave and P. J. Steinhardt, Phys. Rev. Lett. 80 (1998) 1582. 

[7] S.M.Carroll, M.Hoffman and M.Trodden, Phys.Rev. D68 (2003) 023509 (astro-ph/0301273). 

[8] S.Nojiri and S.D. Odintsov, Phys.Lett. B562 (2003) 147-152 (hep-th/0303117). 

[9] A. Kamenshchik, U. Moschella and V. Pasquier, Phys. Lett. B 511 (2001) 265. 
[10] A. Frolov, L. Kofman and A. Starobinsky, Phys.Lett. B545 (2002) 8-16 (hep-th/020418). 
[11] T. Inagaki, X. H. Mcng and T. Murata, hep-ph/0306010. 

[12] E. Elizalde, J. E. Lidsey, S. Nojiri and S. D. Odintsov Phys.Lett. B574 (2003) 1-7 (hep-th/0307177). 
[13] Carroll, S. M, Living Rev.Rel. 4 (2001) 1 (astro-ph/0004075). 



22 



[14] Dvali, G., Gabadadze, G., Porrati, M. 2000, Phys. Lett. B485, 208 (hep-th/0005016); Dvali, G., 
Gruzinov, A., Zaldarriaga, M. 2003, Phys. Rev. D, 68, 24012; Dvali, G., Turner, M. S. 2003 (astro- 
ph/0301510); Dvali, G. Zaldarriaga, M. 2002, Phys. Rev. Lett. , 88, 91303; Freese, K. Lewis, M. 2002, 
Phys. Lett. B540, 1 (astro-ph/0201229). 

[15] A. A. Starobinsky, Phys, Lett. B 91 (1980) 99. 

[16] S. Capozzicllo, S. Carloni and A. Troisi, "Recent Research Developments in Astronomy and Astro- 
physics" -RSP/AA/21-2003 (astro-ph/0303041); S. Capozzicllo, Int.J.Mod.Phys.D 11 (2002) 483. 

[17] Carroll, S. M., Duvvuri, V., Trodden, M. Turner, M. S. 2003, (astro-ph/0306438). 

[18] S.Nojiri and S.D.Odintsov, Phys.Lctt. B576 (2003) 5-11 (hcp-th/0307071); S.Nojiri and S.D.Odintsov, 
Mod. Phys. Lett. A, 19(8) (2004) 627-638 (hep-th/0310045). 

[19] S.Nojiri and S.D.Odintsov, Phys.Rev. D68 (2003) 123512 (hcp-th/0307288). 

[20] D. Barraco, V. H. Hamity, H. Vucctich, Gen. Rcl. Grav. 34 (4) (2002), 533-547. 

[21] X.H.Meng and P.Wang, Class.Quant.Grav. 20 (2003) 4949-4962 (astro-ph/0307354); X.H.Meng and 
P.Wang, Class.Quant.Grav. 21 (2004) 951-960 (astro-ph/0308031); X.H.Meng and P.Wang, (astro- 
ph/0308284). 

[22] D. N. Vollick, Phys.Rev. D68 (2003) 063510 (astro-ph/0306630). 

[23] E.E.Flanagan, Phys.Rev.Lett. 92 (2004) 071101 (astro-ph/0308111); E.E.Flanagan, Class.Quant.Grav. 

21 (2003) 417-426 (gr-qc/0309015). 
[24] R.Dick, (gr-qc/0307052). 

[25] J.P.Duruisscau and R. Kcrncr, Gcn.Rcl.Gra. 15 (1983), 797-807. 

[26] I. Brevik, S. Nojiri, S.D. Odintsov and L. Vanzo, (hep-th/0401073). 

[27] G. Magnano and L.M. Sokolowski, Phys. Rev. D50, 5039 (1994), (gr-qc/9312008). 

[28] T.Chiba, Phys.Lett. B575 (2003) 1-3 (astro-ph/0307338). 

[29] A. D. Dolgov and M. Kawasaki, Phys.Lett. B573 (2003) 1-4 (astro-ph/0307285). 

[30] M.Ferraris, M.Francaviglia and I.Volovich, Nouvo Cim. B108 (1993) 1313 (gr-qc/9303007); M.Ferraris, 

M.Francaviglia and I.Volovich, Class. Quant. Grav. 11 (1994) 1505. 
[31] D. N. Vollick, (gr-qc/0312041). 

[32] A. Borowicc, M. Ferraris, M. Francaviglia and I. Volovich, Gen. Rel. Grav. 26 (7) (1994), 637 
[33] G. Allemandi, A. Borowiec and M. Francaviglia, Acceleration and Inflation in Ricci-square Non-linear 
Gravity, in preparation. 

[34] A. Borowiec, M.Ferraris, M.Francaviglia and I.Volovich, Class. Quant. Grav. 15 (1998) 43-55. 
[35] X.H.Meng and P.Wang, (hep-th/0309062). 

[36] S.Nojiri and S.D.Odintsov, to appear in GRG (hep-th/0308176). 
[37] G.J. Olmo and W. Komp, (gr-qc/0403092) 



23 



